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Abstract 

A dynamical Maxwell system is 

&t = curl/i, ht = —curie in £1 x (0, T) 

e|i =0 = 0, h\t=Q = in fl 

e e = f in 9Hx [0, T] 

where O is a smooth compact oriented 3-dimensional Riemannian 
manifold with boundary, ( • )g is a tangent component of a vector at 
the boundary, e = & (x, t) and h = h* (x, t) are the electric and mag- 
netic components of the solution. With the system one associates a 
response operator R T : f (->■ —v A h' \dUx(o,T)i where v is an outward 
normal to dVl . 

The time-optimal setup of the inverse problem, which is relevant to 
the finiteness of the wave speed propagation, is: given R 2T to recover 
the part Q T := {x G O | dist (x, dQ.) < T} of the manifold. As was 
shown by Belishev, Isakov, Pestov, Sharafutdinov (2000), for small 
enough T the operator R 2T determines fi T uniquely up to isometry. 

Here we prove that uniqueness holds for arbitrary T > and pro- 
vide a procedure that recovers Q T from R 2T . Our approach is a version 
of the boundary control method (Belishev, 1986). 
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Introduction 

0.1 Maxwell system 

Let be a smooth0 compact oriented three-dimensional Riemannian man- 
ifold with the boundary T, g the metric tensor on ft. A dynamical (time- 
domain) Maxwell system is 

e t = curl/i, h t = —curie in int ft x (0, T) (0-1) 

e| t=0 = 0, h\ t=0 = in ft (0.2) 

e e = f onrx[0,T], (0.3) 

where int ft := ft\r, ( • )q is a tangent component of a vector at T, f is 
a boundary control. A solution {e, h} describes an electromagnetic wave 
initiated by the boundary control, e = e*{x,t) and h = h*(x,t) being its 
electric and magnetic components. Since the divergence is an integral of 
motion, the solution satisfies 

div e f ( ■ , t) = div h f ( ■ , t) = in ft (0.4) 

for all t > 0. 

An 'input —> output' correspondence of the system is described by a 
response operator R T : / i— > —v A hs |rx[o,T]> where v is an outward normal at 
the boundary, A is the point-wise vector product. 

A function (eikonal) 

t(x) := dist (x, T) , x G ft (0.5) 

determines the subdomain (near-boundary layer) 

ft T := {x G ft | t{x) < T}. 

By the finiteness of the domain of influence principle for the Maxwell equa- 
tions (shortly: locality principle), the relation 

supp {e f ,h f } C {(x,t) \x e ft, t > t(x)} (0.6) 

holds and provides the exact meaning of that the waves propagate with finite 
speed. 

Everywhere in the paper, 'smooth' means C°°-smooth. 
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By the same principle, the extended problem 



e t = curl h, h t = —curl e 



in D 



XT 



(0.7) 



e = 







h = 



in {(x,t) e D' 
on T x [0, 2T] 



,2T 



t < t(x)} 



(0.8) 
(0.9) 



in a space-time domain 



D 



2T 



{(x,t) | x G intfi T , < t < 2T 



r(x)} 



turns out to be well posed, whereas its solution {e^,h^} is determined by 
the part Q T of the manifold. With the problem fl0.7|) - fl0.9p one associates an 
extended response operator R 2T : / i-> — i/ A h* |rx[o,2T]- 

0.2 Main result 

So, by the locality principle, the operator R 2T is determined by the part fl T 
of the manifold Q. A reasonable question is: to what extent does the 
operator R 2T determine Vt T l As was shown in [2], for small enough 
the operator R 2T determines Q T uniquely up to isometry. Here this result is 
strengthened as follows. 

Theorem 1 For any fixed T > 0, the operator R 2T determines the subdo- 
main Q T up to isometry. 

The proof is constructive: given R 2T we describe a procedure that provides 
a manifold Q T and endows it with a metric tensor g so that (f2 T , g) turns 
out to be isometric to (f2 T , g). 

This result was announced in [B]: as was claimed, it can be obtained 
by straightforward generalization of the approach developed there for the 
acoustical system to the Maxwell system. However, the proof, which we pro- 
pose here, is much simpler, it is based on quite elementary geometric facts 
and the version [8] of the fundamental Holmgren- Joihn- Tat aru theorem on 
uniqueness of continuation of solutions to the Maxwell system across a non- 
characteristic surface. In comparison with the complicated scheme [BJ, which 
uses such devices as the Friedrichs extension and Duhamel integral represen- 
tation, the reconstruction procedure proposed here looks more prospective 
for numerical realization. 

2 namely, for T < dist (c, T), where c is the separation set of w.r.t. T 
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0.3 Comments 

• Setting the goal to determine an unknown manifold from its boundary 
inverse data (here, the response operator R 2T ), we have to keep in mind 
the evident nonuiqueness of such a determination: all isometric man- 
ifolds with the mutual boundary have the same data. Therefore, the 
only relevant understanding of 'to determine' is to construct a mani- 
fold, which possesses the prescribed data [6]. It is what is done in our 
paper: we provide the manifold (f2 T , g), whose response operator R 2T 
coincides with the given R 2T by construction. 

• The literature devoted to inverse problems of electrodynamics is hardly 
observable and we restrict the list of references by the papers dealing 
with the time-optimal setup of the problem given above. Such an opti- 
mality means two things: on the one hand, just by the locality principle, 
no i? 2 ( T_£ ) with e > determines Q T and, on the other hand, to deter- 
mine f2 T we need no R 2{ - T+e " > : it suffices to know R 2T . The longer is the 
time interval of observations at the boundary, the bigger is the part of 
the manifold recovered from the observations. As far as we know, at 
the moment the boundary control (BC-) method is the only approach, 
which provides such a locality of reconstruction %. 

• Dealing with the Maxwell system on a manifold, we use and refer to 
the certain facts and results of [3J, which are proved not for this general 
case but for the Maxwell system in Q C R 3 with the scalar parameters 
e and fi. In all such cases, the generalization is trivial: to get the proof 
one can just put e = fi = 1 and use the intrinsic operations curl, div 
etc relevant to the Riemannian structure instead of the Euclidean one. 

• The authors thank Yu.D.Burago and S.V.Ivanov for very helpful con- 
sultations on geometry. We are grateful to I.V.Kubyshkin for kind help 
in computer graphics. 

• The work is supported by the RFBR grants No 08-01-00511 and NSh- 
4210.2010.1. 



3 In scalar multidimensional inverse problems, the iterative approach by V.G.Romanov 
(see [13], [14]) is also time-optimal. Close results for the Maxwell system see in [15] . 
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1 Geometry 
1.1 Manifold Q 

We deal with a smooth compact Riemannian manifold Q: dim Q = 3, F := 
dfl, d and g are the distance and metric tensor in Q. Also, the notation (Q, g) 
is in use and, if Q is considered as a metric space (regardless its Riemannian 
structure), we write (fi, d). 

Convention In what follows, for a subdomain E C f2, by (E, g) we denote 
this subdomain endowed with the restriction g\-£ of the metric tensor. The 
notation (E, d) means E endowed with the interior distance that is induced 
by the tensor g and measured along the curves lying in E. 
Obviously, (E, g) determines (E, d) and, as is well known, the converse is 
also true: given the distance d on E one can recover the smooth structure 
and the metric tensor g, i.e., determine (E, g). 

The eikonal r( • ) : — >■ R + is defined by (10.51) . For a subset A C Q, we 
denote its metric neighborhood by 

QT [A] := {x e Q\ d(x,A) < r}, r>0; 

in the case of A = T, we write 

Q r ■= Q r [T] = {xeQ \ t(x) < r}. 

Later, in dynamics, the value 

T* :=maxr(-) = inf {r > I fi r = m 

is interpreted as a time needed for waves moving from T with the unit speed 
to fill tt. 

The level sets of the eikonal 

F:={xen \ t(x) = s}, s>0 

are the surfaces equidistant to T. 

For a subset A C Q, by int A we denote the collection of the interior 
points of the set A\T. 
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Figure 1: Cut locus 



1.2 Cut locus 

Recall the definition of a separation set (see, e.g, [10]). Let Z 7 [0,s] be a 
segment of the length s of a geodesic Z 7 emanating from 7 G T orthogonally 
to T; let x(j, s) be its second endpoint. The value r*(7) is said to be a critical 
length if r(x(7, s)) = s for < s < ^(7) (i.e., Z 7 [0, s] minimizes the distance 
between a; (7, s) and T) and r(x(7, s)) < s for s > ^(7) (i.e., Z 7 [0, s] does not 
minimize the distance; see the illustration on Fig 1, where s < 1^(7) < s'). 
Note that the function r*(-) is continuous on T [10J. The point a; (7, t*(7)) 
is a separation point on / 7 . A set of the separation points 



is called a separation set (cut locus) of Q w.r.t. T. 

There is one more way to introduce a cut locus. A point m G Q is said to 
be multiple if it is connected with T through more than one shortest geodesies 
(see Fig 1, where #(7, s) = #(7', s) = m and s = r(m) = ^(7) = ^(7')). 
Denote by Cq the collection of multiple points and define 

c := Co. 

Using (mutatis mutandis) the arguments of [12] , one can check that this 
definition is equivalent to the first one . 
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The cut locus is 'small': the continuity of the function r* easily implies 
vole = 0. 

Note in addition that T s \c is a smooth (may be, disconnected) surface in 
Q. If s < d(c, T) then V s is smooth and diffeomorphic to T. 

1.3 sgc and pattern 

For any x G Q\c, there is a unique point 7(0;) G T nearest to x and a 
pair (7(2), t(x)) is said to be the semigeodesic coordinates (sgc) of x. If 
7 X ,7 2 are the local coordinates in a neighborhood a C T of 7(2), then the 
functions 7 1 (7( • )) > 7 2 (t( ' )) > r ( ' ) constitute a coordinate system in a 'tube' 
{x G I 7(2) G a, < r(x) < n(i(x))}. 
A set 

6 := {( 7 (x), r(ar)) | x G Q\c} = {(7, s) | 7 G T, < s < r,( 7 )} cTx [0,TJ 
is called a pattern of f2, whereas 

:= U^' r *W) cTx [0,T»] 

is its coast that is the graph of the function r*. The sets 

e T : = en {r x [o,t)}, f:=«n{rx[o,T)} 

are referred to as a pattern of Q T and its coast. The patterns are the sub- 
graphs of the continuous functions r*(-) and rj(-) := min {r*(-), T}: see 
Fig. 2@. 

A map i : fi\c — > 0, 

i( x ) : = ( 7 (x),r(x)) 
is a diffeomorphisnjf]. Its inverse i~ l transfers G onto Vl\c by the rule 

« _1 ((7,7-)) = £(7)1-) 

and can be extended to U 9 by the same rule. In the sequel, we deal with 
the extended it maps T U 9 T to Q T continuously and surjectively but 

4 The illustrations on Fig 2 and 3 are taken from [6] 

5 9 is a manifold endowed with the smooth structure of T x [0, T*]. 
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not injectively. As is evident, for two points (7, s), (7', s') G T U 9 T the 
equality 

i - 1 (( 1 ,s)) = t- 1 (( 1 ',s')) (=men T nc) (l.i) 

is valid iff they lie at the coast 6 T , s = s' = r(m) holds, and their mutual 
image m is a multiple point (see Fig. 2). 

1.4 Caps 

Fix 7 G T, s > and a (small) e > 0; let 0^(7) := {7' G T| dfa',^) < e} be 
a portion of the boundary. We say a subdomain 



■■= w M>y)] P| {n-\n-*} = 

= {iG(l d(x,£7 e (7)) < s, s-£<r(x)<s} (1.2) 
to be a cap and note the monotonicity: C Co^' e ' as e < e' . Introduce a set 

lima4' £ := Pi 



w 7 



0<£<S 



and recall that a; (7, s) is defined in sec. 1.2; the following result describes the 
behavior of the caps as e — > 0. 

Proposition 1 T7ie relation 

lim< = { " (7 ' S) if ^ T * (7) (1.3) 

7 \ if s > r*( T ) 1 ; 

holds. 
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Figure 3: Cap 



Proof see in [6] (Lemma 1). 

So, for a given 7 G T and s > 0, the cap either shrinks to the point 
a; (7, s) if (7, s) G U 9, or terminates (disappears for small enough e) if 
(7, s) ^ 8Uft see Fig 3, where the cap is shadowed and ^(7) > s > 
t*(7')) holds. Such a behavior of caps leads to the following evident facts, 
which we'll use for solving the inverse problem. 

Corollary 1 Let T > be fixed. A point (7,5) G T x [0, T) belongs to the 
set G T U 9 T iff for any e > the relation 

^ (1-4) 

holds; in this case, the inequality s < t^(j) is valid. Otherwise, if the family 
of caps terminates, one has (7, s) ^ T U 6 T and, hence, s > t*(j) is valid. 

Corollary 2 Let 7 G T and (7', s) G Q T , so that x' = £(7', s) G Vl T \c. For 
a fixed r < T, the inclusion x' G Q [7] in Q T ( or, equivalently, the inclusion 
(7', s) G i (n, [j)\c\ on Q T ) holds iff the relation 

<f|^ +£ K(T)] ± (1-5) 

i>aM /or any e > 0. 
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Figure 4: Relations (1.5) and (1.6) 

Corollary 3 Le£ i/ie points (7, s) and (7', s) belong to the coast 6 T . The 
equality i -1 ((7, s)) = i^ 1 ((7', s)) i>aM iff for any e > t/ie relation 

u>? f| MY)] ^ (1.6) 

/ioWs (or, equivalently, u^f D f2 s+£ [o" e (7)] 7^ 0j. 

Corollary [T] is just a convenient reformulation of Proposition [TJ whereas Corol- 
laries [2] and |3] easily follow from (jl.3p (see the illustrations on Fig.4a,b, where 
the caps are shadowed). 

1.5 Manifold (Q T , g) 

Here we prepare a fragment of the future procedure that solves the inverse 
problem. The fragment is the following construction. 

Assume that we are given with the pattern T endowed with the tensor 
fi'sgc : = {i~~ l )*9 B As a metric space, the pair (0 T , g sg c) is isometric to the 
space (jl T \c, (f) , d being understood as the interior distance 0. Our goal is 
to make (0 T ,g sg c) into an isometric copy of (f2 T , g). 



6 in other words, <? sgc is the metric tensor g in the semi-geodesic coordinates. 
7 Recall the convention of sec. 1.1! 
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• Step 1 Attach the coast 6 T to the pattern T and extend the metric to 
Q T L>9 T by continuity. On the extended space, introduce the equivalence 

{( 7 , S ) 1 ( 7 ',*')} # {i-\7,8) = r 1 ^,*)} ; (1.7) 

let n T := [<d T U 6 T ] /E be the quotient set, vr : 6 T U 6 T -> fF the 
projection. Recalling the aforesaid about the map i -1 , we see that the 
equivalence class of a (7, s) G T U# T consists of more than one element 
iff i^ 1 ((7, s)) is a multiple point. Also, the map /3 := o n^ 1 : fi T — > 
f2 T is a well defined bijection. So, as result of sewing the proper points 
of the coast, we get a set fl T bijective to fl T . 

• Step 2 Endow T with the quotient topologjj^. As is easy to recognize, 
the bijection f3 is a homeomorphism. Thus, we have a topological space 
Q T homeomorphic to (Q T , g). It remains to endow this space with the 
relevant Riemannian structure. 

• Step 3 Define a set c := ^(0 T ) C f2 T , which is the image of the 
cut locus via bijection Equip the set Q T \c = 7r(6 T ) with the 
metric tensor g : = (ir^)* g agc . As is easy to see, (fi T \c, g) is a manifold 
isometric to (fi T \c, g). 

To extend g to c one can take a point a e c, its neighborhood a; C fi T 
covered by local coordinates it 1 , u 2 , u 3 , find the matrix {gjkiu 1 , u 2 , u 3 )} 3 h 
in u3\c and then extend the matrix elements to lj fl c by continuity. We 
denote the extended tensor by the same symbol g. 

As result, we get a Riemannian manifold (f2 T , g), which is isometric to 
(Q T , g) by construction. 

1.6 Tensor g sgc via distant functions 

Return to the starting point of the procedure (9 T , g sgc ) =>- (Q T , g) and 
explain, where the tensor g sgc will be taken from. Let £l r [y] be a semi-ball 
of the radius r < T with the center 7 e T, z (f2 r [7]\c) C T its image in the 
pattern. 

8 A C f2 T is open iff 7r _1 (y4) is open in Q T U 8 T . 
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Lemma 1 The family of the semi-ball images 

B T := {i (fi r [7]\c) her, < r < T} 
determines the tensor g sgc on T . 

Proof As is evident, to know the family B T is to know the distant functions 

r*(j,T) := d(x(j,T),a) 

for all a G T and (7, r) G T provided the distance in the r.h.s. does not 
exceed T. 

Fix a point (7, r) G T . Take a (small) £ > such that r] e [j, r] : = 
°e(7) x (t — £, r + £:) C 6 T and rd^u,-] < T holds for all a G cr e (j). Let 
7 1 , 7 2 be the local coordinates on (t s (j). The metric tensor g sgc on the pattern 
G T (that is the metric tensor g in sgc) is of the well-known structure 

/ ^s g cii(7 1 ,7 2 ^) ^s g ci2(7 1 ,7 2 ,^) 
9s g c = 9s g c2i(l\l 2 ,r) ^ sg c 22(7 1 ,7 2 , r ) 
\ 1 

=1 we denote the inverse matrix, which is of the same structure. 
Choose three base points Oj G o~ £ {^) (see Fig.5). The equalities 

(^(7, ')) 2 +^ c (7, r) ^( 7 , r) ^(7, r) = 1 i = 1, 2, 3, a, /3 = 1, 2 , 

are just the form of writing the well-known fact: the gradient of any dis- 
tant function on Q is of the norm 1. Consider these equalities as a linear 
system w.r.t. three unknowns g^ c , g^ c = g^ c , g^ c ; let A aia2a3 (7,r) be its 
determinant. As is easy to show, the freedom in the choice of 
quite enough to provide A 0ia2a3 (7,r) 7^ 0. Hence, varying (if necessary) 
the position of the base points, one determines from the system the tensor 
components at the point (7, r). □ 

1.7 Plan 

Here we outline the scheme for solving the inverse problem. 

Assume that the knowledge of the inverse data (operator R 2T ) enables 
one to check the relations (|1.4j) . (11. 5p . and fll.6p . If so, we can realize the 
following construction: 
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Figure 5: Recovering g Sj 



• Select the points (7,5) G T x [0,T), for which (II. 4p holds and thus 
recover the function rji^f), 7 G T, which determines the pattern G T C 
T x [0, T) and its coast 6 T 

• Fix 7 G T, take (7', s) G T , r > available for Corollary [21 and check 
by ( II. 5p whether (7', s) 6 [l]\c) holds. Varying (j'js), recover 
z(n r [ 7 ]\c). Then, varying 7 and r, determine the family B T of semi- 
ball images on T . By Lemma HJ this family determines the metric 
tensor g sgc on 6 T ; thus, we get the manifold (0 T , ^gc). 

• Using (11. 6p as a sewing criterion, introduce the equivalence E. Given 
(0 T , g sgc ) and E, construct the isometric copy (f2 T , g) of the manifold 
(fl T , g) by the procedure of sec 1.5. 

At this point, one can claim that the inverse problem is solved: an iso- 
metric copy of the original manifold is determined from the inverse data. It 
remains to explain, how to check (ll.4l) - (ll.6p via R 2T . 
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2 Dynamics 
2.1 Vector analysis 

Begin with recalling the definitions of operations on vector fields in a 3d- 
manifold (see [16] for detail). As above, g is the metric tensor (2-form on 
vector fields). We assume that fl is oriented and denote by /i the volume 
3-form. 

A scalar product ' ■ ' : {fields} x {fields} — > {functions} is defined by 
a ■ b = g(a, b). 

A vector product A : {fields} x {fields} — > {fields} is defined by g(a A 
b, c) = /i (a, b, c). 

For a field a, its conjugate 1-form a$ is defined by a${b) = g(a,b). Con- 
versely, for an 1-form one defines its conjugate field (jr by (f)(b) = g(b, (jr). 

A gradient V : {functions} — > {fields} is Vu = (dw)", where d is the 
exterior derivative. 

A divergence div : {fields} — > {functions} acts by diva = *d*aj, where 
* is the Hodge operator. 

A curl is defined as a map curl : {fields} — > {fields} , curia = (^dao)'. 

The basic identities are div curl = and curl V = 0. 

Let v be the outward unit normal on T, fir the (induced) surface form 
on T; recall the Green formula 

/ u diva /! = J ua- v — J a ' ^ u • (2-1) 
n r n 



2.2 System a T 

In what follows, we deal not with the complete Maxwell system f l0.ip - (IQ.3p 
but its electric subsystem that is obtained by eliminating the magnetic com- 
ponent 0; we write it in the form 

e tt = -curl curl e in int Q T x (0, T) (2.2) 

e = in{(x,t) E tt T x (0,T) \t < r(x)} (2.3) 

e e = f onrx[0,T], (2.4) 



9 A reason to single out the subsystem is that the components and are not quite 
independent: for times t < T* the magnetic component is determined by electric one [3]. 
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and refer to as the dynamical system a T . Such a form corresponds to the 
locality principle: the solution e* is determined by the part (fl T , g) only 

Let us equip a T with standard control theory attributes: spaces and 
operators. 

Outer space A control / in (12. 4 j) is a time-dependent tangent vector field 
on T. Let T be a space of square summable tangent fields on T; the space of 
controls J-" T := L 2 {[0, T]; T} with the product 

T 

(f,g)TT = J dt J f(j, t) • 0(7, t) nv 

r 

is called outer. It contains an extending family of subspaces 

JT^ . = y e jrT j / | ^ <T _ S = 0} , < C < T 

formed by delayed controls; the parameter £ is regarded as the action time, 
whereas T — £ is the delay. 
Also, the class of controls 

^ : =L 2 ([o,r];rn^(r)) , 

where H?(T) is the Sobolev vector spac^l, and the smooth class 

M T := {/ g r([o,r] ; rnr(r)) |su PP / c (o,t]} 

are in use. Note that by the definition of the second class, each / 6 Ai T 
vanishes near t = 0. The class is a normed space w.r.t. the relevant 
norm, whereas Ai T is dense in J-" T and in J 7 ^. 

For / G A^ T , the problem (I2.2p - fl2.4p has a unique classical (smooth) 
solution e*, which is a time-dependent vector field in f2. Being initially 
defined on A4 T , the map / 1— > e* acts continuously from .7-T to the space 

C ^[0, T], L 2 (^ T ) j (see [2]). Hence, it can be extended to the class J^. In 
what follows we assume that such an extension is done and regard its images 
as generalized solutions of the problem (|2.2p -( |2~4"|) . 

10 as is customary on manifolds, we assume (T) to be endowed with one of the equiv- 
alent relevant Sobolev norms 
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Inner space In control theory, the solution & is referred to as a trajectory 
of the system a T , whereas e*( ■ ,t) is a state at the moment t. By f10.4D . the 
states are divergence-free fields, whereas by (10. 6p one has 

suppe / (-,t) Cff cH T , (2.5) 

— T 

i.e, the trajectory does not leave f2 . Thus, the natural candidate for the 
role of the space of states is 

J T := clos {y e L 2 {^ T ) | supp y C Q T , div y = in Q T } , 

where clos is the closure in the vector fields space L 2 (fi T ) with the product 

(v> v )z2(ffr) = Jyv v 

and divy = is understood in the sense of distributions. The space J T is 
called inner, the waves ( • ,t) are its elements. 

The inner space contains an extending family of subspaces 

:= clos{?/ G J T | supp y C f^} , < ^ < T . 

By the definitions and (I2.5p . one has 

e f (-,t)ej\ 0<t<T. (2.6) 



Control operator An 'input — > state' correspondence in the system a T is 
realized by a control operator W T : J-~ T — > J7" T , DomVT 71 = J-T, 

W T / := e'(-,T), 

which is a closable unbounded operator (see [9]). Hence, it can be extended 
up to a closed operator, what we assume to be done^ and denote the closure 
by the same symbol W T . 

As each closed operator, the control operator can be represented in the 
form of the polar decomposition 

W T = $ T \W T \, (2.7) 

11 However, the precise description of DomlV is not known yet. 
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where \W T \ := [{W T )* W T ]~ 2 is the modulus of W T and $ T is an isometry 
from RanlVT 71 ! onto Raniy T (see, e.g., [7]). 

Response operator An 'input — > output' correspondence is described by 
a response operator R T : J-~ T — > J-" T , Domfl T = Ai T , 

t 

(R T f) (7, t) := A curl J e% s) ds , (7, t) e T x [0, T] , 

'0 

where v is an outward normal, the expression in the right hand side being 
equal to — A/i^|r x [o,r]- This operator is unbounded but closable. The latter 
is a simple consequence of the following fact: it can be shown that R T acts 
continuously from T T to a Sobolev negative space H~ p (T x [0, T]) with a 
big enough p > 0. 

In the mean time, the locality principle enables one to extend the problem 
to the system 

e-tt = —curl curie in D 2T (2.8) 

e = in {(x, t) e D 2T \t <t(x)} (2.9) 

e e = f on T x [0,2T] (2.10) 

(see flO.7p - flO.9p ) and introduce an extended response operator R 2T : J-" 2T — > 
T 2T , BomR 2T = M 2T , 

t 

(R 2T f) (7, t) := 1/(7) A curl y e% s) rfs , (7, t) e T x [0, 2T] . 



This operator is also determined by the part (f2 T , g) of the manifold and, 
hence, can be regarded as an intrinsic object of the system a T 0. 

Connecting form A bilinear form c T : J-" T x — y R, Dom c T = J-T x F + , 



c T [f,g] := (e^(-,T),e^(-,T)) jT = (W T f,W T g) jT 

is called connecting. By closability of iy T , the form c T is also closable [7]. 

The following fact is one of the key points of the BC-method: it is used 
in all of its versions [TJ, [B]. 



12 More about the extended (continued) response operator see in [5]. 
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Proposition 2 The connecting form c T is determined by the response oper- 
ator R 2T . 

Moreover, c T can be expressed through R 2T explicitly as follows. Let S T : 
T 1 - — > T 2T be the operator that extends the controls, as functions of t, 
from [0,T] to [0,2T] by oddness w.r.t. t = T. Introduce the class Ai T '° : = 
{/ EM T \ S T f G M 2T }; it is dense in T% and S T M T '° C Domi? 2T holds. 
The relation 

c T [f,g] = (-l(S T yR 2T S T f,g^ (2.11) 
is valid for / G M T >°,ge 1% (see [3]). 

Corollary 4 The response operator R 2T determines the operator \ W T \. 
Indeed, R 2T determines c T , and one has 

c T [f,g] = (W T f, W T g) jT = (see (E7j)> = (\W T \f, \W T \g) TT , (2.12) 

whereas the positive operator \W T \ is determined by the latter form in the 
r.h.s.. 



2.3 Reachable sets 

Fix an open a C T; a linear set 

*4*M --={W T f\fe M T , supp / C a x (T - £, T] } 

is called reachable (from cr, at the time £). So, the reachable set consists of 
electric waves e*(- ,T) produced by all smooth delayed controls, which act 
from the part of boundary cr, the action time being equal to £. Its closure 
(in J T ) 

is said to be a reachable subspace. Also, we denote shortly W^fT] =: W^. 
For a measurable subset A C fl T , we define a subspace 

J (A) := clos G J" T | supp j/Ci}. 

By the finiteness of the wave speed propagation, each wave belonging to 
^rcach [a] is supported in the metric neighborhood f^[cr]. Hence, the embed- 
ding 

K 6 [o-] C J(tf[a]) (2.13) 
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holds. A structure of the reachable sets and subspaces, as well as the char- 
acter of the embedding (12.131) is the subject of the boundary control theory. 
For our goals, the following fact is of crucial value. 

Proposition 3 Any field h G J~(Q^[a]) QU^[o~] is smooth and satisfies 

curl/i = (2.14) 

in int [a] . 

Proof see in [3], [6] and the remark in the third item of sec 0.3; more 
about properties of the reachable sets see in [4]. Note that the derivation 
of (I2.14p relays upon the fundamental Holmgren- John- Tataru theorem on 
uniqueness of continuation of the solutions to the Maxwell system across a 
non-characteristic surface [8]. 

Corollary 5 If a subdomain uj C intf^[cr] is homeomorphic to an open ball 
in M. 3 , then the relation 

J{u) C U^a] (2.15) 

holds. 

Indeed, by ( I2.14p any h G J(Vfi[a}) U^[a) can be represented locally as 
h\ u = Vp. Hence, for y G J(u) C J {&[a}) one has 

(h,y)jT = J Vp-yyL= (see<$J§) = 0, 

i.e., y G U^[o~] and we arrive at (I2.15p . This result shows that the reachable 
subspace U^[a] is always rich enough. Moreover, in a relevant sense, the 
generic case is that the defect subspace J (Q£[a}) QU^[a] is finite dimensional, 
its dimension being determined by topology of Q^[a] (see, e.g., [3]). 



3 Inverse problem 

3.1 Wave caps 

Here the main technical tool for solving the inverse problem is introduced. 
In parallel with the definition fll.2p . introduce a subspace 



w s f :=w s M 7 )]p|{^ew s - £ } 



(3.1) 
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and call it a wave cap. The correspondence between the wave and 'space' 
caps is as follows. Assume that 7, e, and s are chosen so that the cap u*' £ is 
nonempty. 

• Each y e w^' e is supported in cJ^ 6 . Indeed, by ( 12. 14ft with a = T, the 
orthogonality y-LU s ~ £ implies that the field y is harmonicl^l in £l s ~ £ . 
In the mean time, y is supported in the neighborhood Q be (7)] and, 
whence, vanishes out of this neighborhood. By the uniqueness of the 
harmonic continuation, y has to vanish everywhere in Q s ~ £ and, thus, 
to be supported in the space cap. 

• Fix an x and a (small) ball Q s [x] C intu;*' e . As is easy to see from 
(I2.15p . the subspace J(VL & [x\) is embedded in the r.h.s. of ( 13. ip . Hence, 
the embedding Jj{fl s [x}) C w^ ,£ holds and shows that the wave cap is a 
rich enough nonzero subspace. 

As one can easily conclude from the above-mentioned facts, the space and 
wave caps behave in one and the same way: the equivalence 

{</«}^K £ /{0}} (3.2) 

holds. 

Return to Corollaries 1-3. The 'subdomain <-> subspace' arguments quite 
analogous to the ones that have led to (13. 2p imply the equivalences 

{< f| ^ +£ M7)] * 0} {^f n U r+e We(l)} ^ {0}} (3.3) 

{<n fis+£ w^i ^ } {<n ^ s+£ ^ {o}} • (3.4) 

3.2 Model caps 

Return to the polar decomposition ( 12. 7p . 

Fix an open a C T. In the outer space J-" T , introduce a linear set 

\U\LM ■= {\W\ T f\fEM T , supp/cax (T-&T]} 
13 i.e., curly = and divy = 0; see, e.g., [TB] 
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and a subspace 

\U\*[a) :=clos|W|L ch M = ($WM (3.6) 

(the closure in J-" T ) that we call the model reachable set and subspace re- 
spectively. Subsequently, define a model wave cap 

:= |w| s M 7 )]f|{M s e|w| s - £ } = ($ T )x< £ . (3-7) 

Since the map ($ T )* is an isometry, the subspaces in the right hand sides of 
the equivalences (13 .3p and H3.4[) can be replaced by their model copies, i.e., 
the images through this map. Indeed, for instance, by virtue of the evident 
equality 

H?n i^r^)] = ( $T )* (<' £ n wr+£ [^)0 

the conditions 

<fl ^[^(7)] ^ {0} and H;f H M r+E M7)] ^ {0} 
are equivalent. 

Return once again to Corollaries 1-3. The aforesaid enables one to refor- 
mulate them as follows. 

Corollary 6 Let T > be fixed. A point (7,5) G T x [0,T) belongs to the 
set T U 6 T iff for any e > the relation 

NT 5 ^ {0} (3-8) 

holds; in this case, the inequality s < 1^(7) is waiic?. Otherwise, if the family 
of caps terminates, one has (7, s) ^ G T U 9 T and, hence, s > t*(j) is valid. 

Corollary 7 Let 7 E T and (7', s) G 6 T ; so £/ia£ x' = £(7', s) G fi T \c. For 
a fixed r < T, the inclusion x'gl] [7] in Q T ( or, equivalently, the inclusion 
(7', s) G i (ft [y]\c) on Q T ) holds iff the relation 

f] lUr^)} ^ {0} (3.9) 

holds for any e > 0. 
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Corollary 8 Let the points (7,5) and (7', s) belong to the coast 6 T . The 
relation (7, s) = (7', s) is valid iff for any e > the relation 

f] \K\ s+£ W^)} ^ {0} (3.10) 

holds (or, equivalently, D |W| s+£ [o" e (7)] 7^ {0}^. 

The key fact is that the subspaces in the l.h.s. of (13 .81) — (13 . 1 01) are determined 
by the operator \W T \ and, whence, by the inverse data (operator R 2T ). 

3.3 Reconstruction 

Now, to solve the inverse problem it suffices just to summarize our previous 
considerations, which we present in the form of the following procedure. 
Recall the starting point: we are given with the extended response operator 
R 2T . 

Step 1 Recover the connecting form c T from (12. lip and determine the op- 
erator \W T \ (see Corollary H (T2TT2D ). 

Step 2 Determine the sets |W|j: each [cr], subspaces |W|^[cr], and the model caps 
|w|^' £ by the definitions ( 13.51) . (13.61) . and ( 13 . Tl) respectively for those a, £,7, £, 
which are used in Corollaries EHHJ 

Step 3 Checking the relation ( 13 .8[) . recover the pattern T and its coast 9 T . 
Step 4 By checking ( 13. 9}) . determine the family of ball images B T and 
recover the metric tensor g sgc on the pattern (see Lemma [1]). 

Step 5 By the use of (13.101) . recover the equivalence E. The collection 
Q T ,8 T ,E determines the manifold (f2 T , g), which can be constructed by 
means of the procedure described in sec 1.5. 

As a result, the operator R 2T determines the manifold (f2 T , g), which is 
isometric to the manifold (f2 T , g) by construction. In other words, it deter- 
mines (f2 T , g) up to isometry that proves Theorem 1. Moreover, identifying 
the points by 

n T d r 3 7 = 7T ((7, 0)) g dn T c n T 

we can easily check that the response operator R 2T of the constructed mani- 
fold is identical to R 2T . The latter motivates to refer to (fl T , g) as a canonical 
representative of the class of isometric manifolds possessing the given inverse 
data R 2T . 
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In conclusion, note the following. It is not inconvincible that the above 
described procedure R 2T =>■ (f2 T , ~g) is available for numerical realization. 
The principal problem (and difficulty) is to get good enough simulation of 
the wave caps w^ 6 that is to provide a rich enough set of waves e^(-,T) 
concentrated near the point x(j,s). The sharper is the concentration, the 
better is the reconstruction of the pattern, distant functions r a , etc. Such a 
principle is traced in all versions of the BC-method. 
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